Universal quantum computation with unlabeled qubits 
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I. INTRODUCTION 
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We show that an n-th root of the Walsh-Hadamard transform (obtained from the Hadamard gate 
and a cychc permutation of the qubits), together with two diagonal matrices, namely a local qubit- 
flip (for a fixed but arbitrary qubit) and a non-local phase-flip (for a fixed but arbitrary coefficient), 
can do universal quantum computation on n qubits. A quantum computation, making use of n 
qubits and based on these operations, is then a word of variable length, but whose letters are always 
taken from an alphabet of cardinality three. Therefore, in contrast with other universal sets, no 
choice of qubit lines is needed for the application of the operations described here. A quantum 
i^" , algorithm based on this set can be interpreted as a discrete diffusion of a quantum particle on a de 

' Bruijn graph, corrected on-the-fly by auxiliary modifications of the phases associated to the arcs. 

^ ; PACS numbers: 81P68 
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The study of universality in quantum computation goes back to jj]. In the circuit model, universality has been 
I considered hj a number of papers (see 0, and the references therein) . Probably, the simplest universal set of gates 
consists of the Hadamard gate H together vi/ith the TofFoli gate T "9^. Of course, in order to have universality, we need 
the freedom of applying H and T to arbitrary qubits of the computer: H to any qubit and T to any three qubits. In 
f-->) fact, out of measurement processes, every computational step is induced by a unitary, obtained by tensoring together 
H's, T's and identity matrices. The number of different unitaries obtained in this way is then a function of the number 
. of qubits. 

\o . 

, In this paper, we define a universal set of unitary quantum operations depending on the computational space of the 
machine. This means that the unitaries change whenever the total number of qubits changes. Since the operations 
Oh! act globally on all qubits, they do not require the choice qubit lines at each computational step. This is in contrast to 
^Jj ' standard finite universal sets {e.g., H and T) where the gates reman fixed but may be applied to varying choices of 
Ch qubit lines. The set defined here is composed by an n-th root of the Walsh-Hadamard transform (constructed from 
the Hadamard gate and a cyclic permutation of the qubits), and two diagonal matrices: a (local) qubit- flip (for a 
fixed but arbitrary qubit) and a (non-local) phase-flip (for a fixed but arbitrary coefficient). To prove that the defined 
set is universal on n qubits, we reduce it to H and T, acting on any qubit and any three qubits. The characteristic 
property of the set is that the choice of the qubits to which apply H and T is not reflected into the structures of 
the unitaries. This observation would like to be a justification to the title of the paper. Of course, the set described 
; ^ > here is inconvenient from the physical impementation point of view, because it requires nonlocal interaction between 
' qubits (which is physically expensive). However, it may be useful to remark that the set provides a mathematical 
framework, in which doing universal quantum computation is conctructing words whose letters are unitaries (two 
of which commute) taken from an alphabet of cardinality three. So, the outcome of a computation depends on the 
length of the word and the order of the letters. This has some flavour that reminds of quantum finite state automata 
and other sequential machines (see e.g. |3)I3|)- 



The remainder of this paper is organized as follows. In Section 2, we give some preliminary definitions. In Section 
3, we formally state and prove the main result. This is done by proving universality with a reduction to {i?, T}. It 
is not difficult to verify that the n-th root of the Walsh-Hadamard transform respects the topology of the de Bruijn 
graph. Namely, the ij-th entry of this unitary is nonzero if and only if there is a directed edge from the vertex labeled 
i to the vertex labeled j, in the de Bruijn graph on 2" vertices. In Section 4, in virtue of this observation, we point 
out that any quantum algorithm can be seen as the discrete diffusion of a quantum particle on a de Bruijn graph, 
corrected on-the-fly by a qubit-flip and a phase-flip, both fixed but arbitrary. This reminds of the context of discrete 
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quantum walks [31 or the processes studied in p. A natural open question would be to prove that the n-th root of 
the Walsh-Hadamard transform and a phase-flip (for a fixed but arbitrary coefficient) form a universal set. 



II. DEFINITIONS 



In this section we introduce some preliminary definitions. 

Definition of Vn- We denote by Vn a square matrix of dimension 2" such that G {0,±-^} and with exactly 

the following nonzero entries: 



[Vn 
[Vn] 



[Vn] 



1,2"-1 



[Vn 



_ J_ 

" V2' 

J2,2"-1 + 1 = -75, 
[Ki]3,2 = [Ki]3,2"-i+2 = [Ki]4,2 = y| , 
[K]4,2"-i+2 = ^75' 



(1) 



[Ki]2'»-1,2"- 
[K]2",2" = 



— [K]2"-l,2" — [Ki]2",2"-1 — 7f , 

'72- 



The matrix Vn is real-orthogonal and it is an n-th root of _ff„ := iJ®", where H is the 1-qubit hadamard gate: 



(2) 



For example, 



Vo 



1 

V2 



/ 1 1 

I -1 
10 1 

Vol 0-1 



(3) 



and V2 = H . An alternative and more direct definition of Vn can be given as follows. Let Sn be the full symmetric 
group on the set {1, 2, n}. We denote permutations of length n as ordered sets. For example, the elements of S'3 
are (1, 2, 3), (1, 3, 2), (3, 2, 1), (2, 3, 1), (2, 1, 3), and (3, 1, 2). With an abuse of our notation, their regular permutation 
representations are denoted in the same way. The matrix Vn is defined as 



Vn :-P-(i/®/®"-i), 

where 

P=(l,3,...,2"-1,2,4,...,2"). 
This permutation is nothing but a cyclic-shift of the qubits: 

P : \aia2...an) — > |a„_iai...a„_2)- 



(4) 
(5) 
(6) 



This explains why — Hn {Hn '■= -ff®"). It is easy to verify that Eq. ^and Eq. 2]define the same matrices. Notice 
that, when n = 2, the permutation P is the Swap-gate: 



(1,3,2,4) 



/ 1 
10 
10 

Vo 1 



(7) 



Definition of Pn{k). Let us denote by P„(fc) the 2" x 2" matrix defined by 



if i ^ j; 
[Pn{k)]i^] { -1 if i = fc with k e {1,...,2"}; (8) 

1 otherwise. 
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The matrix Pn{k) is the Pauh operator Z = 



1 

-1 



on the rt-th qubit controlled by all other qubits 1, 2, rt — 1 



and then conjugated by the permutation that interchanges dimensions k and 2". For example, 

^2(3): 



10 
10 
0-10 
1 



= /2 © xzx, 



(9) 



where X is the Pauli operator X ~ 



1 

1 



Definition of Fn{k). Let us denote by F„(fc) the 2" x 2" matrix defined by 

r if ^ ^ j; 

[Fn{k)]i^j := < —1 if the fc-th qubit is 1; 

I 1 otherwise. 

The matrix Fn{k) is the Pauli operator Z acting on the fc-th qubit and the identity on all other qubits. For example, 



(10) 



F2{1) 



10 
0-100 
10 
-1 



/(g) Z and F2(2) = 



/I 








1 











-1 





\0 





-1 



^Z(g>I. 



(11) 



III. MAIN RESULT 

In this section we prove the following theorem. 
Theorem 1 The set 

B = {Vn,Pn{i),F^ij)}, 

for fixed but arbitrary i e {1, ...,2"} and j € {1, is universal for quantum computation on n qubits. 



(12) 



By this theorem, any quantum computation in the circuit model can be seen as a word whose letters are taken from 
the alphabet B. 

The Toffoli gate T is the 3-qubit gate defined as T : (a, 6, c) — > (a, 6, a6 © c), where © denotes addition modulo 2 
and a,b,c G {0, 1}. The matrices defined by the Toffoli gate are then a special kind of transposition. If we label the 
states of the computational basis in lexicographic order, the matrices defined by the Toffoli gate on 3 qubits are 



/10000000\ 
01000000 
00100000 
00010000 
00001000 
00000100 
00000001 

VOOOOOOlO/ 



/10000000\ 
01000000 
00100000 
00000001 
00001000 
00000100 
00000010 

Voooioooo/ 



and 



/10000000\ 
01000000 
00100000 
00010000 
00001000 
00000001 
00000010 

VOOOOOlOO/ 



(13) 



The set B' — {H,T} is universal for quantum computation We prove Theorem ^ by showing that the set 

B is equivalent to the set B' , once the number of qubits has been fixed. In order to verify this equivalence, we show 
that any expression made by the tensor product of H, T, and the 2x2 identity matrix I2 corresponds to a sequence 
of the three elements of B defined in the statement of the above theorem. We will consider the set 



Dn := {Pnii) : 1 < * < 2"}. 
Lemma 2 The set A ^ VnU Dn is universal for quantum computation on n qubits. 



(14) 
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Proof. By the definition of V,i (Eq. 0] above) , we have 
where 

P-I = (1, 2"-i + 1, 2, 2"-i + 2, 3, 



$n-l 



(15) 



S2") = (l,3,...,2"-l,2,4,...,2")-i. 



This means that if we can construct the matrix P ^ then we can also construct the Hadamard gate. Firstly, observe 
that 



y-l ^ y2n-l ^ yT 



Since 



The matrix can be then obtained directly from y„. The nonzero entries in the last two rows of Vn are 



1 

V2- 



[14]2",2" 

The nonzero entries in the last two columns of Qn = P,i(2") • Vj[ are 



Since H-^{\ \ 



[Qn]2"-1,2"-1 — [Qri]2",2"-1 — [Pn(2") • Kf]2"-1.2" 
[Qri]2",2" = ^"Tf 

X . it follows that 



1 



(16) 



VnQn = KPn(2")l^J = /2"-2 © X = (1, 2, 2" - 2, 2", 2" - 1), 

which is indeed the Toffoli gate. Now, given that 

^2" = ((l,2,...,2"-2,2",2"-l),(2,3,...,2",l)), 

if we can construct the cyclic permutation (2, 3, 2", 1) then we will have the full symmetric group 82^- If we can 
construct 52" then we will get and all the permutations that we need for applying H and T to arbitrary qubits. 
The permutation (2, 3, 2", 1) can be constructed with the following procedure (most probably not optimal): 

1. Let Ml = HnFn{n)Hn = By applying elements from I?„, tranform the bottom-left block of Mi into 
the diagonal matrix dia{pi, ...,p2n-3, — P2"-3+i, ~P2n-2), where pi = (1, —1). Let M2 be the matrix obtained 
in this way. 

2. Let M2 — M'2Vn- By applying elements from D„, transform the top-right block of M2 into the identity 



matrix and the bottom- left block into the diagonal matrix dia{pi, ■■■,P2" 



'P2" 



3. Repeating the second step n ~ 2 times one gets the permutation matrix (2,3, 
checked with any computer algebra system.) 

An example with three qubits may help to clarify the procedure: 

i/3 (1, 2, 3, 4, 5, 6, 7, 8) i?3 = (5, 6, 7, 8, 1, 2, 3, 4) , 

Fg-i (5, 6, 7, 8, 1,2, 3, 4) 1/3 = (3,4,5,6,7,8,1,2), 

Vf' (3,4,5,6,7,8,1,2)1/3 = (2,3,4,5,6,7,8,1). 



4 + 1, -p2"-3j- 

.,2",1). (This can be easily 



The notation is easily explained: 



/ 



(5,6,7,8,1,2,3,4) 



10 
0-100 
0-10 

V 1 



10 
10 
10 
1 
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By the above constructions, we have the following fact: the matrix group G = (Kt, D„) has the regular permutation 
representation of 82^ as a subgroup. Then T, obtained with Eq. 1161 can be applied to any three qubits; H, obtained 
from Eq. 1151 can be applied to to any qubit. Since B' = {H,T} is universal, the lemma follows. ■ 

Lemma 3 The matrix group G = {VmPn{i),Fn{j)) , for fixed hut arbitrary i G {1, -..,2"} and j G {1, ■■■,n}, contains 
the set Dn- 

Proof. By the definitions, 

= P-^ 

= I®---®I®X®I---®I, 

where X is at the fc-th position of the tensor product. For example, V3'^P3(3)V3~"'^ — I ® I ® X . Let H be the matrix 
group generated by the matrices I®---®I®X,I®---®I®X®I,X®I®---®I. The group H is isomorphic to Z2 
(indeed, the above matrices are the permutation representations of the standard generators of ). Since, for every 
*j j G {Ij 2"}, there is an element of H sending P„(z) to Pn{j), we can construct that is the set of all 2" x 2" 
diagonal matrices with entries 1 and —1. We considered Fn(n), but we could also take P„(j), for any j S {1, ---jn}, 
without loss of generality. ■ 

Theorem ^ is a consequence of Lemma |21 together with Lemma |21 

IV. A RELATION WITH DE BRUIJN GRAPHS 

Let E be an alphabet of cardinality d and let S* be the set of all words of length n over E. The d-ary n- dimensional 
de Bruijn graph is a directed graph denoted by B{d,n) and defined as follows j^j: the set of vertices is E^; there is 
an arc from i to j if and only if the last n — 1 letters of i are the same as the first n — 1 letters of j. The graph 
B{2, n) is called (directed) binary de Bruijn graph. Notice that P(2, n) has exactly two loops: one loop is at the vertex 
. . . 0; the other one at the vertex 1 ... 1. These graphs have important applications in cryptography and distributed 
computing. In particular, they provide some of the best-known topologies for communication networks. For example, 
the Galileo space probe of NASA used a network based on a de Bruijn graph to implement a signal decoder |8|. Let 
M„ be the adjacency matrix of P(2, n). The rows and the columns of this matrix can be ordered in such a way that 

[M„]i,i = [M„] 1^2-1 + 1 = [Mnh,! = [Mn]2,2"-i + l 1> 
[Mn]3,2 = [M„]3_2--i+2 = [Mn]4,2 = [M„]4_2"-i+2 = 1, 

[M„]2„_i,i2-.-l = [Af„]2"-1,2" = [M„]2„^l2„_i = [M„]2.,2. = 1. 

Tanner |llj pointed out that the matrix Vn is obtained from the matrix Af„ by negating the entries 
[M„]2,2"-i+i' [-^^n]4,2"-i+2i [Afn]2",2" and rcscaling M„ by It is easy to see that a simple random walk on 
B{d, n) converges very quickly to uniformity, in fact it is perfectly mixed after n steps. Simple means that at each 
vertex the walker chooses to cross an incident edge by tossing a fair die with d faces. We associate the vertices of 
B(2,n) with the elements of the computational basis |0) = |0...0), |1) = |0...01), |2"-i) = |l...l). For every \i) and 
there is a diagonal matrix Q g D„ such that HnQHn\i) = V^QV^\i) = We may interpret this process as a 
discrete quantum walk on P(2, n) induced by Vn and corrected on-the-fly by an appropriate diagonal unitary with ±1 
entries. If the walk is not corrected then V^"|i) = \i) for every i, given that is symmetric. The walk on B{2,n) 
is perfectly mixed at the n-th step, but it can be driven with probability 1 to any vertex in exactly 2n steps. For 
example, the following is an algorithm that takes the state |0) to the state |2" — 1) in 2n steps: 

V;m = |+)«"; 
Pi(i)®"|+)®" = IV'); 
K"|V) = |2"-1). 
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(A curiosity: the positions of the minus sign in the state correspond to the numbers with an odd number of I's 
in their binary expansion (A007413 10])). A generahzation for any \i) and |j) is straightforward. 

Aknowledgments. I would hke to thank the anonymous referees. Their comments helped me to improve the 
presentation of this paper. 



[1] D. Aharonov, Simple Proof that Toffoli and Hadamard are Quantum Universal, quant-ph/0301040. 

[2] A. Ambainis, Quantum walks and their algorithmic applications, quant-ph/0403120. 

[3] D. Deutsch, Quantum computational networks, Proc. Roy. Soc. Land. A 425 73-90, 1989. 

[4] S. Gudder, Quantum automata: an overview, Internat. J. Theoret. Phys. 38 (1999), no. 9, 2261-2282. 

[5] M.-C. Heydemann, Cayley graphs and interconnection networks. Graph symmetry (Montreal, PQ, 1996), 167-224, NATO 

Adv. Sci. Inst. Ser. C Math. Phys. Sci., 497, Kluwer Aead. Puhl, Dordrecht, 1997. 
[6] J. Kosik, Scattering Quantum Random Walk, Optics and Spectroscopy 99, 224-226 (2005). 

[7] C. Moore, J. P. Crutchfield, Quantum automata and quantum grammars, Theoret. Comput. Sci. 237 (2000), no. 1-2, 
275-306. 

[8] B. Mukherjee, Optical Communication Networks, Series on Computer Communications, McGraw-Hill, New York, 1997. 

[9] Y. Shi, Both Toffoli and controUed-Not need little help to do universal quantum computation, quant-ph/0205115. 
[10] N. J. A. Sloane, (2005), The On-line Encyclopedia of Integer Sequences, www.research.att.com/~njas/sequences/. 
[11] G. Tanner, Spectral statistics for unitary transfer matrices of binary graphs, J. Phys. A 33 (2000), no. 18, 3567-3585. 



